The role of boson-fermion correlations in the resonance theory of superfluids 
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Correlations between a composite boson and a fermion pair are considered in the context of the 
crossover theory of fermionic to bosonic superfluidity. It is shown that such correlations are the 
minimal ingredients needed in a many-body theory to generate the right boson-boson scattering 
length in the Bose-Einstein limit of the crossover. 
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In the quantum theory of many-body systems, super- 
fluidity is connected with a broken symmetry which is 
not anticipated from the underlying microscopic Hamil- 
tonian. There are two distinct types of superfluidity. A 
first type occurs in fermionic systems in the presence 
of an effectively attractive fermion-fermion interaction 
(as in superconducting metals). The symmetry breaking 
there is associated with the development of mean-field 
pair-correlations with coherence length much larger than 
the interparticle spacing. A second type of superfluid- 
ity occurs in bosonic systems (such as atomic conden- 
sates) where the symmetry breaking is associated with 
the condensation of preformed bosons with spatial size 
much smaller than the interparticle spacing. For many 
years it has been conjectured that these two types of 
superfluidity are in fact alternative limits of a single uni- 
versal phenomenon, and that one could continuously pass 
from fermionic to bosonic superfluidity by properly vary- 
ing the fermion interaction parameters. 

This conjecture has attracted a great deal of attention. 
It has been shown that the Bardeen-Cooper-Schrieffer 
(BCS) theory Q reduces to Bose-Einstein condensation 
(BEC) of bosonic dimers as the pair correlation length 
becomes small compared with the interparticle spac- 
ing 0. For a dilute Fermi gas, in which the mean in- 
terparticle spacing is large compared to the extension of 
the two-body potential, the simplest way to obtain the 
crossover is via a scattering resonance. This corresponds 
to bringing the highest-lying bound state from just above 
to just below threshold. In a dilute gas the atom-atom 
interaction can be parametrized by the only nonvanish- 
ing contribution to the two-body T-matrix at low en- 
ergy, namely the fermionic s-wave scattering length ap 
(with T = 4:TTh 2 aF/'m and m is the atomic mass). The 
crossover from fermionic to bosonic superfluidity thus 
appears through the tuning of ap from negative val- 
ues, through infinity, to positive values. A series of 
improvements to the basic theory to better account for 
the crossover region have been developed and extensively 
studied during the last decades 

It has recently been pointed out 0, 0| that although 
these theories model the passage from fermionic to 
bosonic superfluidity, they have one major flaw in that 
they fail to reproduce the correct equation of state of the 



system in the bosonic limit. This is due to the fact that 
they predict the wrong value for the boson-boson scat- 
tering length, cib — 2aF- A solution of the four-fermion 
Schrodinger equation for contact scattering in vacuum, 
which in principle can be determined exactly, has recently 
been found and provided the value as — O.Gap. 

Progress in experiments with quantum atomic gases 
has lead to experimentally accessible systems which can 
directly address the crossover region using Fcshbach scat- 
tering resonances Bose-Einstein condensation of two- 
fermion dimers has been observed and the experimental 
study of the crossover is actively being pursued. The 
first data collected on the Bose-Einstein condensed cloud 
were consistent with the 0.6ap J]. A consequence is 
that crossover theories that predict in the BEC limit the 
value 2a f for the boson-boson scattering length will not 
yield the correct results for all observables dependent on 
interactions, e.g. the collective modes, the vortex core 
structure, the internal energy, and even the macroscopic 
spatial extent of a confined cloud. 

The aim of this paper is to present a theory of the 
crossover which correctly reproduces the Gross-Pitaevskii 
and Bogoliubov theory of the non-ideal Bose gas with the 
right boson-boson interaction. As we shall see this the- 
ory will recover on the fermionic side the BCS theory 
including the Gorkov corrections due to density and spin 
fluctuations @- We begin by pointing out the essential 
ingredients which are missing in standard crossover the- 
ories and which have to be included to achieve this goal. 
Our approach was motivated by Ref. || where the few- 
body problem of four fermions interacting in vacuum was 
considered. The essential property leading to the correct 
boson-boson T-matrix was identified within the contact 
scattering formalism as the three-particle correlation be- 
tween a composite boson and a fermion pair. This fol- 
lowed from the observation that in the limit in which the 
distance r\ between any two of the fermions becomes 
small, the four-fermion wavefunction must factorize as 



V(r u r 2 ,R) -» /(r 2 , J R)(l/4 7 rr 1 - 1/4™*-), 



(1) 



where r 2 is the distance between the two other fermions 
and Rj y/2 the distance between the centers of mass of the 
pairs. The function /(r 2 ,-R) — an effective three-body 
wavefunction of a composite boson and two fermions — 
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contains all information needed to determine the boson- 
boson scattering length. The reason why standard 
crossover theories fail to give the correct T-matrix in the 
BEC limit is then revealed, because they are typically 
based on the BCS assumption that only two-fermion cor- 
relations are important, and only those are withheld in 
the correlation hierarchy. 

In this paper we show that a self-consistent many-body 
theory of the crossover can be constructed which includes 
three-particle correlations, and correctly reduces to the 
known results in the appropriate limits. This modifica- 
tion not only brings about quantitative corrections, but 
also introduces a qualitative change in the crossover pic- 
ture. A major consequence is that the superfluid order 
parameter in this crossover theory is not given by the 
simple integration of a BCS-type pairing field. 

A system of fermions in two internal states (f and J.) 
close to a Feshbach resonance can be described using the 
many-body Hamiltonian 0, Il0| : 



where e k = h 2 k 2 /2m is the free fermion dispersion rela- 
tion, <7fc is the matrix element relating two free fermions 
in the open channel to the closed channel bound state 
near threshold, and v is the bare detuning of the bound 
state. The operators a^l annihilate (create) open chan- 
nel fermions with momentum k and spin <r, while 6^ an- 
nihilate (create) closed channel bosons. We ignore here 
the direct interaction of the fermionic degrees of freedom 
which would otherwise give rise to an asymptotic back- 
ground value for the scattering length, since that plays no 
role at detunings sufficiently close to the resonance value. 
This can nevertheless be reintroduced in a straightfor- 
ward manner as has been previously discussed [llj . 

In resonance superfluidity theory the information con- 
tained in the wavefunction /(r2,R) for a boson and a 
fermion pair is in principle conveyed by the correlation 
function 



( h -1%/2- 



-k[ a q/2+k 



(3) 



Notice, however, that the two quantities do not simply 
coincide because the one in Eq. © describes the corre- 
lations between bare fermionic and bosonic states, while 
/(r 2 , R) gives the correlations between states dressed by 
the interactions (and thus physically observable). The 
connection between the two quantities is not obvious and 
must be stated explicitly. 

The relation can be well understood by considering the 
two-fermion correlations with zero center of mass mo- 
mentum, which in a vacuum coincides with the Fourier 



transform of the relative wavefunction. The bare corre- 
lation function is simply \&-kl®kt/- This can however 
never be an eigenstate of Eq. j2j ■ An eigenstate can be 
constructed by considering the following linear combina- 
tion of (fi_/yflfei-) and (b Q ) (dressed state): 



*fc = (a- fci a fcT ) + V 



9k 



2e k -E 



(4) 



with V denoting the Cauchy Principal Value and E a 
solution of 



at 



2e k -E 



(5) 



The nature of the solution depends on the presence or 
absence of a bound state indicated by the sign of the 
renormalized detuning v. This is defined as 



(6) 



and is physically related to the magnetic field shift from 
the Feshbach resonance ^lj. The case of v < cor- 
responds to the BEC side of the resonance. There a 
bosonic dimer bound state exists and the solution of 
Eq. JSJ coincides with the bound state energy E = 
-h 2 /ma 2 F [H [H Q. For V > there is no bound 
state and the solution is E = v. Evolving Eq. Q for 
two particles in a vacuum under the Hamiltonian Eq. © 
gives the Schrodinger equation 



ih 



dt 



2e fe * fe + E U k ,k'^k' 



with real separable potential 



k.k' 



= V 



9k9k' 
2e k -E' 



(7) 



(8) 



What remains is to show that this potential generates the 
correct scattering length at all detunings v. To this end 
we must obtain the T-matrix by solving the Lippmann- 
Schwinger equation as derived from Eq. J7J 0] 



T, 



k.k' 



u, 



k.k' 



E 



2e k i — 2e q + iS' 



H 



(9) 



in the limit of zero scattering energy. Substituting Eq. (JSJ) 
into Eq. @ and performing the integration, we recover 
the desired result for the zero-energy T-matrix 



T 



9 2 o 



(10) 



Eq (|10|l provides the correct behavior of the tuning of the 
scattering length around resonance [Tlj . with the usual 
definition T = inh ap /m, and confirms that the po- 
tential in Eq. JSJ leads to the correct effective fermion 
interaction properties. 
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The key point is that Eq. Q is also the time- 
dependent Schrodinger equation for the relative wave- 
function \6t-kl®kt) °f tw0 P ar ticles in a vacuum evolving 
under the single-channel Hamiltonian 

Single = e fe"L a fe<7 

ho 

+ Uk > k ' A l/2+k1; & l/2-kl 6: q/2-k>i a q /2+k>-\ C 11 ) 

qkk' 

Introducing the dressed wavefunction Eq. is therefore 
equivalent to eliminating the bosonic degree of freedom 
from the theory by introducing effective purely fermionic 
quantities. We wish to point out that there is an alterna- 
tive method of eliminating the bosonic degrees of freedom 
based on the application of the Hubbard-Stratonovich 
transformation to the functional integral expression for 
the partition function This is distinct from our ap- 
proach, because it leads to a formal result in which the 
effective fermion-fermion interaction potential contains 
bare quantities, namely v. In contrast the potential we 
give here is a real effective potential which contains only 
renormalized quantities, E or v. 

The implication of the above considerations goes be- 
yond what we have presented so far, because the dressing 
procedure can be extended to treat correlation functions 
of any number of particles. It is now possible to make 
the direct link between the bare and the dressed correla- 
tion functions for a composite boson and a fermion pair, 
Eq. © and /(r 2 ,i?). This can only be done assuming 
the contact scattering form for the boson-fermion vertex, 
i.e. that g k = g independent of k. Since this assumption 
leads to an ultraviolet divergent theory, we need to intro- 
duce a momentum cutoff K to the wavevector sums ^lj . 
This implies for instance that the renormalized detuning 
in Eq. © is given by v — v — mKg 2 /2it 2 Ti 2 . Recur- 
sively applying the same procedure used for the dressing 
of the two-fermion wavefunction, we are led to construct 
an ansatz for the four-fermion wavefunction: 

*fcl,fc 2 .q = ( a - q /2-k 1 i a -q/2+k 1 ] a q/2-k 2 i a q/2+k 2 '\) 
+/3fel ( b -q a q/2~k 2 i a q/2+k 2 l) 

+Pk 2 {b q a_ q/2 _ ki[ a_ q/2+ki] ) 

-P{q-k 1 -k 2 )/2{b hl _ k2 a_ q /2-k l i a q /2+k 2 ]) 

-P(q+k l +k 2 )/2(b k2 _ ki a q/2 _ k2l a_ q/2+ki1 ) 
+PkJjk 2 {b- q b q ) 

-(3(q-k 1 -k 2 )/2f3(q+k 1 +k 2 )/2(b k2 _ kl & fcl _fc 2 > (12) 

where (3 k = 'P{g/(2e k — E)}. It can be shown that for 
four-particles in vacuum the evolution of this ansatz un- 
der the resonance Hamiltonian Eq. (J2J is equivalent to 
the evolution of the four-particle wavefunction 

( a -q/2-k 1 i a -q/2+k 1 '\ a q/2-k 2 i a q/2+k 2 ^ ( 13 ) 



under the single channel Hamiltonian Eq. i|ll|l . with the 
potential U k k > given in Eq. JSJ. Since we have proved 
that this interaction leads to the correct two-body T- 
matrix, this result is crucial because it is exactly the same 
equation which was solved in Ref. [j| to obtain the correct 
value of the boson-boson scattering length. It should 
also be clear that the procedure we have outlined can be 
systematically extended to incorporate more and more 
particles. 

We are now in a position to make a definitive statement 
about the complete set of correlation functions in the res- 
onance theory needed to reconstruct the information rep- 
resented by /(r2,i?). Having two fermions approaching 
each other is equivalent to summing over one of the two 
relative momenta, e.g. fei, in Eq. or equivalently in 
Eq. 1|12|) . In this way we formulate the complete set: 

(b- q b q ) , 

> (14) 

E fel (a_ g/2 _ feli a_ g/2+felT v 2 _ fe2i a g/2+fe2T >. 

The desired many-body theory must therefore include 
these correlation functions in order to generate the cor- 
rect equation of state in the bosonic limit. We can now 
proceed to construct an approximate many-body the ory 
using the standard methods of cumulant expansion |l7j | 
so that we keep explicitly the correlation functions in 
Eq. (|14|) and all the correlation functions of the same 
order. The key point is that the resulting theory closes 
within a set of functions which depend on a maximum of 
two vector arguments (e.g. the functions given in Eq. Ijl4(l 
depend at most on q and fe 2 )- At a computational level, 
this is an essential point, because a two vector field in a 
translationally invariant system contains three nontriv- 
ial degrees of freedom, and a three dimensional theory is 
tractable. This is the minimal complexity because it is 
consistent with the intrinsic dimensionality of f(r2,R) 
which encapsulates the few-body scattering calculation 
that this theory must reduce to in the limit of zero den- 
sity. We also point out that this is an enormous sim- 
plification over the result of a direct application of the 
cumulant expansion to keep all correlation functions of 
order four fermions and below, i.e. including those of the 
form given in Eq. I |13(l . Such a function depends on three 
vectors, and even in a uniform system, where rotational 
symmetries can be exploited, contains at least six non- 
trivial degrees of freedom. 

It must be emphasized that the possibility of reducing 
the number of degrees of freedom in the whole crossover 
region is peculiar to the resonance Hamiltonian. This is 
because in order to close the set of equations in the many- 
body theory with three-particle correlation functions, it 
is crucial to apply the contact scattering approximation. 
In the resonance theory, Eq. this implies g k — > g in- 
dependent of k. However, the resulting effective potential 
Eq. (JSJ) has a residual fc-dependence which makes the po- 
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tential and thus the theory well-defined even at v = 0. 
However, in the single channel formulation, Eq. this 
point causes problems as the scattering length ap — ► oo 
and the two-body T-matrix diverges. For this region 
there is no obvious way of incorporating correctly cor- 
relations of more than two particles in the many-body 
theory in a pseudopotential approximation. One can of 
course utilize a nonlocal potential and keep the full four- 
particle functions of the form in Eq. in that case. 

We now consider explicitly the consequence of the 
above formalism in the bosonic and fermionic superflu- 
idity limits. The complete theory of the dilute nonideal 
Bose gas in the bosonic limit emerges naturally since it 
depends only on the eigensolution of the two-body prob- 
lem and the interactions of the resulting dimers. Both of 
these elements we have discussed in depth. Thus, when 
the usual approximations are made, the Gross-Pitaevskii 
equation, the Bogoliubov theory, and the Popov theory, 
will emerge. In the fermionic limit, the set of correlation 
functions that we are keeping leads to the possibility of 
constructing a pair wavefunction that correctly includes 
the full many-body effects. In particular one can prove 
that the following ansatz: 

*fe = (a-fel^t) _ mAr/ 2 +fe((Ml+feT a fet) 
q 

-(b- q a_ H al q _ H )), (15) 

evolves using Eq. J2J as (oi-kl^kl) ev °l ves under the 
Hamiltonian in Eq. 1|11JI : 

ih = 2e fe (a_ fci a fcT ) 

~ H U q/2+k,k> {i & \+k] & k] & q/2-k' 1%/2+k'^) 



qh' 



*'T>)' ( 16 ) 



kl a ~q/2-k'i a -q/2+ 



Eq. I|16fl shows that Eq. (|15fl contains the full pair wave- 
function in the medium and not just the BCS approxi- 
mate factorized form. Notice that an analogous ansatz 
and correspondence can be formulated for the density 
correlation function (o4|d fc |). Together these imply that 
we not only recover the complete BCS equations but also 
the effects beyond mean-field such as the Gorkov correc- 
tions to the superfluid gap. 

In conclusion, a consistent theoretical framework has 
been presented for the crossover of superfluidity from the 
fermionic to the bosonic type. The inclusion of the cor- 
relation function between a boson and a fermion pair 
was required to reproduce the correct bosonic equation of 
state. This has lead us to develop a fundamental change 
to the picture of the crossover physics. 
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